function [btTp,StTp]=ksmooth(H,F,shat,sig,sfor,sback,omega,omegaback)

J = sig*(H*F)'/(H*omega*H');
btTp = shat + J*H*(sback-sfor);
StTp = sig + J*H*(omegaback - omega)*H'*J';
